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Arterial spin labeling characterization of cerebral perfusion during 
normal maturation from late childhood into adulthood: normal 
'reference range' values and their use in clinical studies 

Patrick W Hales\ Jamie M Kawadler\ Sarah E Aylett^'^ Fenella J Kirkham^ and Christopher A Clark^ 

The human brain changes structurally and functionally during adolescence, with associated alterations in cerebral perfusion. We 
performed dynamic arterial spin labeling (ASL) magnetic resonance imaging in healthy subjects between 8 and 32 years of age, to 
investigate changes in cerebral hemodynamics during normal development. In addition, an inversion recovery sequence allowed 
quantification of changes in longitudinal relaxation time (Ti) and equilibrium longitudinal magnetization (Mq). We present mean 
and reference ranges for normal values of Ti, Mq, cerebral blood flow (CBF), bolus arrival time, and bolus duration in cortical gray 
matter, to provide a tool for identifying age-matched perfusion abnormalities in this age range in clinical studies. Cerebral blood 
flow and Ti relaxation times were negatively correlated with age, without gender or hemisphere differences. The same was true for 
Mo anteriorly, but posteriorly, males but not females showed a significant decline in Mq with increasing age. Two examples of the 
clinical utility of these data in identifying age-matched perfusion abnormalities, in Sturge-Weber syndrome and sickle cell anemia, 
are illustrated. 
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INTRODUCTION 

The transition from childhood, through adolescence and into 
adulthood, marks a period of structural and functional change in 
the normal development of the human brain. The cerebral cortex 
undergoes its most dynamic structural change in adolescence,^ 
and previous studies have indicated that this is accompanied by a 
further rapid change in cerebral perfusion (subsequent to the 
dynamic changes that occur in the first year of life), which occurs 
between the ages of 10 to 20 years.^ 

The nonionizing nature of magnetic resonance imaging (MRI), 
coupled with its superior spatial resolution, makes it an attractive 
technique for studying the developing physiology of the brain. 
Although metabolism is not measured directly with MRI, cerebral 
perfusion can be quantified, which under physiologic steady-state 
conditions is coupled to the level of cerebral oxygen and glucose 
consumption,^ and as such mirrors the brain's metabolic demand 
and neuronal function. 

Traditionally, MRI-based perfusion techniques have involved 
serial imaging after the injection of a paramagnetic contrast agent, 
to dynamically track the passage of a labeled bolus through the 
vasculature. However, the injection of an exogenous contrast 
agent limits the applicability of these techniques in longitudinal 
studies, and precludes their use in patients with renal failure, due 
to the associated risk of nephrogenic systemic fibrosis."^ Its 
application in children has also been limited due to reduced 
patient comfort, and technical difficulties in administering the 
intravenous contrast agent. 



Arterial spin labeling (ASL) is an emerging alternative MRI 
technique for fully noninvasive quantification of cerebral blood 
flow (CBF).^'^ Although the signal-to-noise ratio (SNR) is inherently 
lower than measurements acquired using injected paramagnetic 
contrast agents, ASL has the advantage that it uses an 
endogenous tracer, by magnetically labeling water in the arterial 
blood supply. As such, it is a completely noninvasive and 
nonionizing technique, allowing for safe repeated measurements 
and increased patient comfort. Furthermore, ASL provides 
measurements of CBF in fully quantitative physiological units 
(such as mL blood/1 00 g tissue per minute), and as such is well 
suited for use in longitudinal or multicenter studies. These 
combined benefits have seen ASL move from the field of 
research into routine clinical practice in recent years.^'^ 



Arterial Spin Labeling Perfusion Quantification 
Perfusion quantification using ASL is performed by inverting the 
longitudinal magnetization of the arterial blood flowing into the 
tissue, waiting for a given 'inflow time' (Tl), then acquiring an 
image (known as the 'label' acquisition). The same process is 
repeated without labeling the inflowing arterial blood (the 
'control' acquisition), after which the perfusion information is 
contained in the subtraction between the two acquisitions (dM, in 
which dM = control — label signal intensity). 

A source of potential systematic error when estimating CBF 
using ASL is the delay between the application of the label and 
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the arrival of the labeled bolus of blood at the imaging voxel 
(known as the bolus arrival time, BAT). Spatial variations in BAT will 
result in a variable amount of tagged blood having been delivered 
at the time of the acquisition (Tl), and failure to account for this 
can lead to significant errors in CBF quantification. 

A second source of error is the 'through-flow' effect, which 
occurs when the acquisition is made at an early Tl, leading to an 
overestimation in CBF, as the dM signal includes a significant 
contribution from intravascular tagged blood which is destined to 
perfuse more distal tissue. Both the above problems are linked to 
the fundamental issue that the time required for labeled blood to 
travel from the tagging region to the capillary bed is both spatially 
variable and generally comparable to the relaxation time of 
blood, ensuring that multiple dynamic processes are occurring at 
the time of acquisition (delivery, exchange, clearance by Ti 
relaxation, and flow).^ 

A simplified kinetic model was introduced by Buxton et aP^ to 
account for some of these processes, in which the time- 
dependent difference signal (for pulsed-ASL) is given by: 



dM(f) 
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Here, a = inversion efficiency, Mqb = equilibrium longitudinal 
magnetization of blood, Tib = longitudinal relaxation time of 
blood, Tl = longitudinal relaxation time of tissue, A = equilibrium 
tissue/blood partition coefficient of water, and T= bolus duration. 

If measurements are made at multiple TIs, then the dM time 
series can be fit to equation (1), with BAT, CBF, and T as fitted 
parameters. However, due to SNR and scan time constraints, the 
majority of clinical studies currently acquire ASL data at a single Tl. 
As shown in equation (1), the amount of labeled blood entering a 
voxel will increase linearly between BAT<f<BAT + r, and as such, 
the optimal Tl should be >BAT + r, to eliminate the effect of 
transit time on CBF quantification. A challenge when designing an 
ASL protocol is that the delay required for this is not generally 
known a priori. To help counteract this, it has been shown that 
'Q2TIPS'" saturation pulses can be applied after a given delay (Th) 
after the labeling, to 'cut-off the trailing end of bolus, ensuring a 
well-defined bolus width leaves the labeling region. However, the 
value of Twill also be influenced by bolus dispersion between the 
labeling region and the imaging voxel, which varies as a function 
of vessel size, length, and tortuosity. As such, it would be 
advantageous to know what starting values and local variations in 
BAT and T to expect when designing a single-TI ASL protocol. 

A further point to note is that, for accurate quantification of CBF, 
Mob, Tib, and Ti must be known. A series of inversion recovery 
acquisitions (with varying TIs) can be performed to measure Ti 
and Mo on a voxelwise basis. From this, Mqb can be determined, 
using Mob = Mo/A^^ (note calculating this on a voxelwise basis 
accounts for radio frequency coil Bi inhomogeneities). Due to 
difficulties in locating a voxel filled purely with blood, literature 
values of Tib ^tre generally used, which in studies involving 
children and young adults, should be adjusted for age and 
sex.^^'^^ However, in many previous studies, the 'q{ty term in 
equation (1) has been assumed to be unity, and the effects of both 
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the venous clearance and the exchange of labeled blood water 
into the surrounding tissue have been ignored. As the clearance of 
tagged water by flow is much slower than the rate of delivery, the 
former assumption is generally valid in the healthy human brain. 
However, it has been shown that exchange between capillaries 
and tissue appears to occur at TIs > 1,000 ms after the application 
of the tag.^^ As single-TI measurements are generally made at 
Tl> 1,000 ms, the magnitude of the error introduced by the latter 
assumption will depend on the difference between Tib and the 
local Tl value. In normal gray matter (GM), this is generally small; 
however, in white matter and under certain pathologies, this error 
can be significant. Also, Ti in cortical GM is known to vary with 
age,^^ and as such failure to account for this will lead to a variable, 
age-dependent error in CBF quantification. 

Due to the rapid alterations in CBF that occur during the 
transition from late childhood into adulthood,^ it is important that 
age-related changes in both the dynamic ASL signal and the Ti 
and Mo values used to calculate CBF are fully understood in this 
age group. Although a number of studies have been published in 
which ASL has been used to measure cerebral perfusion as a 
function of age, the majority of these only report CBF values, 
generally because the data were acquired with either a single or a 
few TIs,^'^^"^^ and as such the changes in the dynamic ASL signal 
(and associated changes in BAT and T) have not been extensively 
studied, particularly in children and young adults. In this study, we 
investigated how Ti, Mq, and hemodynamic properties of the 
dynamic ASL signal vary as a function of age, in healthy subjects 
between 8 and 32 years of age. The lower limit of this age group 
was chosen to ensure all subjects could be scanned without the 
use of sedation, as different kinds of sedation may modify CBF in 
different ways,^° and previous studies have shown significant 
increases in CBF in halogen-sedated compared with awake 
children.^ We provide 'reference range' limits (see Materials and 
Methods) for normal values in both raw ASL data and fitted 
parameters for nonsedated subjects, to be used both as a tool in 
future clinical studies and when designing an optimized single-TI 
ASL protocol. We also show two clinical examples of how these 
reference range values can be used to identify age-matched 
perfusion abnormalities: in one patient with Sturge-Weber 
syndrome and one patient with sickle cell anemia (SCA). 

MATERIALS AND METHODS 

Subjects 

Approval to conduct this study was obtained from the Research Ethics 
Committee at Great Ormond Street Hospital. All healthy subjects provided 
written informed consent, in accordance with our institutional ethical review 
board. A total of 30 healthy subjects (7 males and 9 females below 18 years, 
7 males and 7 females above 18 years) were studied, aged between 8 and 
32 years (mean 18.5 years). In addition, ASL was performed in conjunction 
with the standard clinical imaging protocols in use in our institution in two 
patients. The first was a 17-year-old Sturge-Weber patient with a right pial 
angioma, who was experiencing symptomatic focal epilepsy, migraine 
headache, learning difficulties, and intracranial hypertension. The second 
was a 12-year-old patient with homozygous SCA with no evidence of 
neurologic abnormality seen on T2-weighted MRI. Blood taken in this 
patient 6 weeks before MRI revealed a hemoglobin level of 8.6g/dL (normal 
range: 1 1 .5 to 1 5.5 g/dL) and hematocrit of 23% (normal range: 35% to 45%). 
Peripheral blood oxygen saturation was taken as 96% on the day of MRI. 

Magnetic Resonance Imaging 

All experiments were performed on a 1.5-T Siemens Magnetom Avanto 
scanner (Siemens, ErIangen, Germany), equipped with 40mT/m gradients 
and a 12 channel head receive coil. Arterial spin labeling data were 
acquired using a flow-sensitive alternating inversion recovery pulsed-ASL 
sequence, with 3D single shot GRASE data acquisition,^^ with the following 
imaging parameters: repetition time = 3.0 seconds, echo time = 31.6 ms, 8 
averages, field of view = 230 mm, matrix size = 64x64, 20 contiguous 
slices with 5 mm thickness. Measurements were made at six TIs, ranging 
from 0.2 to 2.2 seconds in 0.4-second intervals, with a total scan time of 
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4.8 minutes. Background suppression of static tissue was used, as 
described by Gunther et al,^^ with x 10 scanner gain, to maximize the 
signal from the inflowing blood. A Ti-mapping sequence was then 
performed in each subject, comprising a series of inversion recovery 
acquisitions (Tl = 0.2, 0.6, 1.4, and 2.4 seconds) with identical readout, field 
of view and resolution as the ASL sequence, but without the use of 
background suppression. 

Postprocessing 

All data analysis was performed using Matlab (MathWorks Inc., Natick, MA, 
USA), and all model fitting was performed using an iterative Nelder-Mead 
nonlinear least squares algorithm. 

First, the GM in the ASL and Trmapping images was segmented into 
major vascular territories, using the flow territory templates established by 
Tatu et al}^ Following these templates, regions of interest (ROIs) defining 
the left and right cortical vascular territories of the anterior cerebral artery 
(ACA), middle cerebral artery (MCA), and posterior cerebral artery (PCA) 
were drawn on a standard MNI (Montreal Neurological Institute) template. 
These were then registered onto each subject's ASL/Ti -mapping images, 
using an affine 1 2 parameter model provided by the FLIRT algorithm in FSL 
(FMRIB's Software Library, Oxford University, www.fmrib.ox.ac.uk/fsl). These 
ROIs were manually corrected in each subject, where necessary, to correct 
for registration errors caused by the point spread function inherent in the 
3D GRASE acquisition. An example from a representative subject is shown 
in Figure 1. 

After image segmentation, regional values of T-i and Mq were calculated 
by fitting the mean signal intensity from the 7"i-nnapping acquisition in 
each ROI to the following equation: 

5(0 =Mo[l -2e"^ + e"?] (2) 

with It and Mq as fitted parameters. A mean value of Mqb was then 
calculated in each ROI, using Mob = Mo/A, with A = 0.9.^^ After this, the 
mean ASL dM(f) signal from each ROI was fit to equation (1), using the 
calculated values of It and Mqb, and assuming a = 1. Age- and gender- 
corrected values of Tib (at 1.5T) were used for each subject, based on the 
data presented in Piechnik et al?^ This provided fitted values of BAT, CBF, 
and T in each ROI. In addition, both the peak value of the fitted dM curve 
(d^max) and time at which the peak occurred (Tl^ax) were recorded. 

Reference Range Regions 

Using the method described by Royston,^^ data from all healthy subjects 
were used to calculate both correlations and reference ranges for normal 
values of each fitted parameter as a function of age. This technique 
provides an age-dependent mean value and a 'normal range' for a given 
parameter, accounting for the fact that both the mean and the standard 
deviation of the raw data may vary with age. Data were log transformed 




Figure 1. Exannple of arterial spin labeling (ASL) innage segnnenta- 
tion, based on vascular flow territories. An axial ASL difference 
innage (dM) is shown, at inflow tinne (Tl) = 1.8 seconds, with a 
sennitransparent overlay illustrating six vascular territories. L, R = left, 
right; ACA, MCA, PCA = anterior-, nniddle-, and posterior-cerebral 
arterial flow territories, respectively. 



where necessary to ensure normality of the residuals, and a backwards 
stepwise regression method was used, starting with a second-order 
polynomial, then reducing the degree to a first-order polynomial (i.e., linear 
fit), provided a significant worsening of the fit was not observed (details in 
ref. 24). All reference ranges were calculated with 95% confidence intervals. 

Statistical Analysis 

All comparisons between correlations were performed using an analysis of 
covariance (ANCOVA) test, followed by a Tukey-Kramer post hoc analysis 
(significance was defined as P<0.05). As part of the post hoc analysis, 
population marginal means (PMMs) were compared, which represent the 
predicted value of the dependent variable for the mean value of the 
independent variable (age in our case). Differences in PMMs are listed in 
the Results section in cases where significance was reached (P<0.05). 



RESULTS 

In all ROIs, none of the fitted parameters showed a significant 
difference in their correlation with age as a function hemisphere 
(i.e., left- and right-hand sides were equivalent). Therefore, all 
regional values represent the mean for the left-/right-hand side of 
the brain. A significant difference between males and females was 
only seen in one parameter (Mq in the PCA territory, see J^/Mq 
Values section). In all other cases, a single correlation/reference 
range is provided, for male and female subjects combined. 

Ti/Mo Values 

The variation in and Mq as a function of age is shown in 
Figure 2. Ti showed a significant negative correlation with age 
(P< 0.001) in all regions. The PMM was significantly lower for Ti 
values in the PCA territory (1 .1 7 seconds), compared with both the 
MCA and ACA territories (1.23 seconds and 1.24 seconds, respec- 
tively). Mo also showed a negative correlation with age in the ACA 
and MCA vascular territories. The vascular territory of the PCA was 
one case in which a significant difference between males and 
females was observed. Here, males showed a highly significant 
negative correlation between Mq and age (P< 0.001); however, no 
significant correlation was observed in females (P = 0.11). In this 
case, data are shown separately for male and female subjects 
(Figure 2F). 

dMr^ax/Tlmax Values 

The variation in dMmax and TImax as a function of age is shown in 
Figure 3. The magnitude of the peak ASL signal (dMmax) showed a 
strong significant negative correlation with age in all vascular 
territories (P< 0.001). In these regions, the mean dM^ax value at 8 
years of age was over double the mean value at the age of 30. 
There were no significant differences between the PMM dM^ax 
values in all regions (P = 0.11). The Tl at which the peak occurred 
(TImax) showed no correlation with age in all vascular territories. 
The PMM was significantly higher on the PCA territory (1.74 sec- 
onds) than that found in both the ACA (1.57 seconds) and MCA 
(1.51 seconds) territories. 

Bolus Arrival Time, CBF, and 7 Values 

The variation in BAT, CBF, and T as a function of age is shown in 
Figure 4. Fitted values of BAT showed no correlation with age in all 
vascular territories, and the PMM value in the PCA (0.53 seconds) 
was significantly higher than that found in the ACA (0.40 seconds) 
and MCA (0.41 seconds) territories. Cerebral blood flow showed a 
negative correlation with age in all regions, which was particularly 
strong in the ACA territory (P = 0.0005) compared with the MCA 
(P = 0.03) and PCA {P = 0.02) territories. The PMM value was higher 
in the MCA territory (63.4mL/100g per minute), compared with 
both the ACA (56.4mL/100g per minute) and PCA (54.2mL/100g 
per minute) territories. A significant positive correlation between 
age and 7 was observed in the PCA (P = 0.0005) territory; this fell 
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Figure 2. Variation in fitted values of Ti and Mq as a function of age. Data are shown for the cortical vascular territories of the anterior (left 
colunnn), nniddle (center colunnn), and posterior (right colunnn) cerebral arteries. Data points represent nnean values for left- and right-hand 
sides of the brain. Data are shown for nnales and fennales connbined, except in panel F, in which a significant difference in the correlation 
between Mq vs. age was observed in nnales connpared with fennales; here, individual correlations are shown (black circles = nnale, gray 
triangles = fennale). Ti values are shown in the top row (A-C) and Mq values are shown in the bottonn row (D-F). Reference ranges are 
indicated in gray (for a 95% confidence level), with the fitted curve illustrated in black. The order of the polynonnial fit is indicated by 'poly' (i.e., 
poly= 1 is linear and poly = 2 is quadratic), and t indicates that the data were log-transfornned before fitting. The significance of the fit is 
indicated by: ^**P< 0.001, **P<0.01, *P<0.05. 



just short of significance in the ACA (P = 0.07) territory, and no 
correlation was observed in the MCA territory (P = 0.4). The PMM 
value of T was significantly higher in the PCA territory 
(1.32 seconds), connpared with both the ACA (1.22 seconds) and 
MCA (1.1 4 seconds) territories. 

The raw values for the nnean curves, and upper and lower 
reference range linnits, for all fitted paranneters shown above, are 
given in the Supplennentary Material. 

Examples of the dynamic ASL signal, calculated using 
equation (1) and incorporating the mean, age-matched values of 
BAT, CBF, r, Ti, and Mq, for representative 10- and 30-year olds, are 
shown in Figure 5. Between 10 and 30 years of age, the MCA 
territory showed the largest decrease in the magnitude in the ASL 
signal, however, the overall shape of the dM curves was 
approximately similar for the two ages. Conversely, the PCA 
territory showed the smallest difference in the magnitude of the 
ASL signal, but the dM curve at 30 years was much 'wider' and 
'flatter' compared with the curve at 10 years (this was also evident, 
to a lesser extent, in the ACA territory). This results in a greater 
fitted value of T. 



Clinical Example in a Sturge-Weber Syndrome Patient 
Figure 6 shows example r2-weighted and postgadolinium Tr 
weighted images for a patient with Sturge-Weber syndrome, along 
with the fitted ASL and Ti -mapping parameters which fell outside 
the reference range regions in this patient. The Tj- and Ti -weighted 
images revealed cerebral atrophy and the corresponding leptome- 
ningeal enhancement on the patient's right-hand side. Elevated 



cortical Ti values (above the reference range) were found through- 
out the right cerebral hemisphere (RCH), particularly in the vascular 
territories of the PCA and the ACA. Reduced CBF was observed 
throughout the brain in the RCH, and was also below the reference 
range in the ACA territory in the left cerebral hemisphere. Elevated 
values of T were observed bilaterally in the ACA territory. 

Clinical Example in a Neurologically Asymptomatic Sickle Cell 
Anemia Patient 

Figure 7 shows example T2- and Ti -weighted images for the SCA 
patient included in this study, along with the fitted ASL and Ti- 
mapping parameters that fell outside the reference range regions 
in this patient. Two independent neuroradiologists concurred this 
patient has no evidence of cerebral infarction on r2-weighted MRI. 
Reduced cortical Ti values were observed throughout the right 
hemisphere, as well as on the left-hand side in the MCA territory. 
Very high values of CBF were observed bilaterally in all regions 
(142, 145, and 130mL/100g per minute in the vascular territories 
of ACA, MCA, and PCA, respectively, mean of left- and right-hand 
sides). These values were approximately double the average value 
in healthy age-matched subjects. Reduced values of T were 
observed bilaterally, in vascular territories of ACA and MCA. 

DISCUSSION 

Healthy Subjects 

Our CBF measurements in normal GM agree well with a number 
of previous ASL studies. For instance, in a cohort with a mean 
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Figure 3. Variation in the nnagnitude of the peak arterial spin labeling (ASL) signal (dMmax/' A-C), and the inflow tinne (Tl) at which the peak 
signal occurs (Tl^ax/ D-F), as a function of age. All synnbols are identical to those used in Figure 2. 
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Figure 4. Variation in the fitted values of bolus arrival tinne (BAT, A-C), cerebral blood flow (CBF; D-F), and bolus duration (7; G-l) as a function 
of age. All synnbols are identical to those used in Figure 2. 



age of 1 1 ± 3 years (nnean ± s.d.), Jain et al^^ reported a nnean well with our mean value at 11 years of 62±4mL/100g per 
CBF value of 65 ± 10 mL/IOOg per nninute, calculated using minute (mean of all GM regions ± s.d.). At the older end of our 
measured, subject-specific blood Tl values, which agrees very cohort, our mean value of 50 ±6 mL/IOOg per minute at 32 years 



Journal of Cerebral Blood Flow & Metabolism (2014), 776-784 



© 2014 ISCBFM 



ASL MRI characterization of cerebral perfusion 
PW Hales et al 



781 




Figure 5. Exannples of the typical dynannic arterial spin labeling (ASL) signal (d/W) in representative healthy 10- and 30-year olds (red and blue 
curves, respectively). The dM curves were calculated using the nnean, age-nnatched values of bolus arrival tinne (BAT), cerebral blood flow 
(CBF), T, 7"i, and Mq shown in Figures 2 to 4. Data are shown for the vascular territories of the anterior cerebral artery (ACA; A), nniddle cerebral 
artery (MCA; B), and posterior cerebral artery (PCA, C). 




' left » right 

middle pDSterior 
1 4r 




10 20 30 Ifl 20 30 10 20 30 



E 

§ 50 



10 20 30 10 20 30 10 20 30 



E 2 




10 20 30 10 20 30 

Age(vrs) 

Figure 6. Exannple standard clinical innages and arterial spin labeling (ASL) perfusion paranneters acquired in a 17-year-old Sturge-Weber 
patient. Axial r2-weighted and postgadoliniunn Ti-weighted innages are shown in panels A and B, respectively. The ASL and Ti-mapping 
paranneters that fell outside the reference ranges in this patient are shown in panels C-E. The reference range acquired in healthy subjects is 
shown in gray in all plots, for vascular territories of the anterior- (left-hand colunnn), nniddle- (center colunnn), and posterior- (right-hand 
colunnn) cerebral arteries. Ti values are shown on the top row, cerebral blood flow (CBF) on the center row, and T on the bottonn row. 
Measured values in this patient are overlaid, and are shown for each hennisphere (left = green data points and right = red data points). 



agrees well with the study by Wang et al,^^ who nneasured 
a mean value of 58 ± 1 2 nnL/1 00 g per nninute in an adult group 
with nnean age of 32 years. It should be noted however that the 
range of CBF values in nornnal GM reported in the literature 
varies considerably, as nneasured values will be influenced by 
factors such as the chosen TI in single-TI studies, and often fixed, 
adult values of blood TI are used for data processing, which will 
lead to an overestinnation of CBF in younger subjects 
(particularly girls, due to differences in hennatocrit levels 
during development). Alternatively, positron emission tomo- 
graphy is considered by many to be the current standard in 
human CBF measurements. Considering a subset of the data 
presented in Leenders et aF^ (including only their subjects 
between 22 and 32 years), using ^^O positron emission 
tomography, GM CBF was measured as 59±17mL/100g per 
minute at 27 ±4 years, which compares well with our value of 
53±5mL/100g per minute at the same age. Finally, Chiron 
et al^^ using ^^^Xe SPECT, measured mean GM CBF values of 



62 ± 1 0 mL/1 00 g per minute in the 6 to 1 9 years age range (no 
mean age is given), and 51 ± 8 mL/1 00 g per minute in the 1 8 to 
29 years range, which also agrees very well with our values. 

The strong negative correlation between dMmax and age shows 
how certain aspects of pediatric physiology provide a natural 
solution to overcome the limitation of low SNR inherent in the ASL 
technique. The increased ASL signal at a younger age is due in 
part to increased blood flow in children; our CBF measurements 
show that between the ages of 8 to 32 years CBF showed a 
significant negative correlation with age in all regions. This 
contradicts a previous computed tomography study by Winter- 
mark et al,^^ in which regional CBF levels were shown to peak 
between 2 and 4 years of age, and then stabilize at adult levels at 
~7 years of age. Alternatively, using ^^^Xe single-photon emission 
computed tomography, Chiron et aP^ measured increasing levels 
of CBF from birth until 5 or 6 years of age, followed by a steady 
decrease, reaching typical adult levels (defined as those found in 
19- to 29-year olds) between 15 and 20 years (as data were 
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Figure 7. Example standard clinical images and arterial spin labeling (ASL) perfusion parameters acquired in a neurologically asymptomatic 
12-year-old sickle cell anemia patient. Axial r2-weighted and Ti-weighted images are shown in panels A and B, respectively. The ASL and 7i- 
mapping parameters that fell outside the reference ranges in this patient are shown in panels C-E. All symbols are identical to those used in 
Figure 6. 



averaged for subjects over 19 years in this study, later age-related 
changes in CBF may have been lost). Both Wintermark's and 
Chiron's studies involved radiation and so included children with 
an indication for assessment but in whom parenchymal imaging 
was normal, whereas ASL may be used in completely normal 
children. A previous ASL-based study by Taki et aP^ showed a 
range in GM CBF trajectories between 5 and 18 years of age, 
varying between negative linear to cubic relationships, depending 
on the location. The age at which CBF peak values occurred also 
varied by region (from <5 years in the occipital lobe to ~ 1 1 years 
in the frontal lobe). Finally, our results show the best agreement 
with a previous ASL study by Biagi et al,^ who showed elevated, 
constant CBF levels between 4 and 10 years, which then 
decreased rapidly during adolescence, until a plateau was 
reached around 25 to 30 years of age. 

Our data confirm that Ti relaxation times in GM tissue show a 
strong negative correlation with age between 8 and 32 years. This 
further improves the SNR of the ASL signal in younger subjects, as 
following exchange at the capillaries, the longitudinal magnetiza- 
tion of the labeled blood water will decay according to of the 
surrounding tissue. As this is longer in younger subjects, the 
effective half-life of the tracer is prolonged. The increased Tib 
values at younger ages^^ contribute further to this increase in SNR. 
A similar decrease in Mq with age was seen in most regions. The 
age-related reductions in both and Mq are likely to be a result of 
both increased water content of the brain in children compared 
with adults^^ and the onset of cortical thinning that occurs during 
adolescence.^^ It is interesting that the evolution of Mq in the PCA 
territory is the only region in which a significant difference 
between males and females was seen, with males showing a 
significant decline in Mq with increasing age, whereas females 
showed no significant correlation. The data presented in ref. 29 
showed that cortical thinning occurs at different rates in males 
and females during adolescence, depending on the location. This 
could be the underlying cause for the different rates of change in 
these parameters between males and females in the PCA territory; 
however, this remains speculative, as high-resolution anatomic 
imaging was not performed as part of this study, so GM cortical 
thickness could not be measured. 



Previous transcranial Doppler ultrasound studies have shown a 
significant decline in time-averaged flow velocities (TAV) in ACA, 
MCA, and PCA, between children aged 10 to 17.5 years, and adults 
aged 20 to 63 years.^° Between children aged 1.5 to 9.9 years and 
children aged 10 to 17.5 years, no change in TAV was observed in 
ACA and PCA, but a decline was seen in the MCA. In general, TAV 
was highest in the MCA, followed by the ACA, and lowest in the 
PCA. However, in this study we observed no change in BAT as a 
function of age in all vascular territories. As a higher TAV should 
lead to shorter values of BAT, this would suggest that TAV is 
constant in all vascular territories between 8 and 32 years of age. 
However, these results should be interpreted with caution, as the 
temporal resolution of the ASL measurements made in this study 
is likely to be too coarse to identify small-scale changes in the 
underlying BAT. Furthermore, although some voxels within each 
ROI will include signal from large arteries branching from ACA, 
MCA and PCA, the majority of voxels will be distal to these feeding 
arteries, and hence the ASL parameters predominantly reflect 
blood flow in arterioles and capillaries, which will be less sensitive 
to changes in the TAV upstream in the feeding artery. However, 
the fact that that the PMM value of BAT was significantly higher in 
the PCA vascular territory, compared with both the MCA and ACA 
territories, does support previous findings that TAV is lowest in the 
PCA. Furthermore, our finding that CBF was significantly higher in 
the MCA territory compared with both the ACA and PCA territories 
agrees with previous evidence that TAV is highest in the MCA. 

Around the age of 8 years, all three regions (ACA, MCA, and PCA 
territories) showed similar values of T (1.1 seconds in our case, 
although absolute values of T will vary depending on the 
acquisition protocol and the use of QUIPSS-II pulses). However, 
the relative change in 7 as a function of age progressed differently 
in each region, remaining constant in the MCA territory, while 
showing a highly significant increase with age in the PCA territory, 
and a trend (P = 0.07) for increasing values with age in the ACA 
territory. As shown in Figure 5, the increased value of 7 was due to 
a 'wider/flatter' dM curve in these regions in older subjects, 
particularly in the PCA territory. The increase in T suggests an 
increasing level of bolus dispersion in these territories with 
increasing age (note, an increase in the overall volume of tagged 
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blood with age, due to slower blood velocities in the carotid 
arteries, could also be a cause of longer values of 7, as QUIPSS-II 
pulses were not used in this study. However, if this was the case, 
then a change in 7 would have been observed in all regions). The 
general kinetic model used in this study makes the common 
assumption that the labeled bolus can be modeled using 'plug- 
flow' in the vasculature, in which there is no dispersion. However, 
it has been shown that this is not representative of fluid dynamics 
in an ASL acquisition, and the labeled bolus is both pulsatile and 
has a time varying, nonuniform cross-sectional flow profile.^^'^^ 
A number of models have been suggested to correct for this, for 
instance by introducing a Gaussian distribution of tag arrival times 
at the imaging voxel, and assuming parabolic flow in major 
arteries.^^ As the temporal width of the dM curve seems to 
increase with age, particularly in the PCA territory, it would appear 
that the labeled blood bolus is becoming more dispersed with 
increasing age, and these models would therefore become 
particularly relevant in this vascular territory. Furthermore, this 
natural increase in bolus dispersion with increasing age should be 
accounted for when assessing changes in T in the presence of 
pathologies in future clinical studies. 

Clinical Examples 

The aim of this study was to provide reference range values to 
describe the hemodynamic properties of the ASL signal in healthy 
children and young adults, which can be used as a tool in future 
studies to investigate pathologies that lead to abnormal cerebral 
perfusion. Sturge-Weber syndrome is a good example of one such 
pathology, characterized by intracranial venous dysplasia, a 
cutaneous capillary angioma, and ocular abnormalities due to 
anomalous embryonic development.^^ An abnormal, tortuous pial 
vasculature is present, associated with venous stasis, leading to 
perfusion defects, which have been observed both focally in 
regions of leptomeningeal disease and in more distant GM in 
ipsilateral and contralateral regions.^^ Clinically, there is commonly 
epilepsy and contralateral hemiparesis in addition to high 
intraocular pressure (glaucoma). The Sturge-Weber syndrome 
patient included in this study had a right pial angioma, with 
evidence of leptomeningial enhancement in the right PCA 
territory (identified on Ti -weighted images, postcontrast), and 
atrophy in the RCH. Ti values were consistently higher in the RCH 
in all vascular territories, which may be due to the effects of 
edema. Cerebral blood flow was below the reference range in all 
regions of the RCH, and also in the left ACA territory. This 
illustrates how hypoperfusion can be observed both locally and 
remote to regions of leptomeningial disease, and also the 
importance of accounting for age when considering CBF values. 
If this patient's CBF values had been compared with healthy 30- 
year-old adults, then the ACA and PCA values would appear within 
the normal range, however, the values are low for a 17-year-old 
patient. The decreased values of CBF are likely to be due to 
impaired venous clearance in these regions. The elevated values 
of T bilaterally in the ACA territory may be representative of 
regions of increased tortuosity in the vasculature. It is interesting 
that this occurs in regions contralateral to leptomeningial 
enhancement; previous studies have shown perfusion abnormal- 
ities before the onset of 'structural' abnormalities indicated by 
conventional imaging,^^ which may be the case in this patient. 
Follow-up imaging would be required to confirm this. 

Sickle cell anemia is another pathology in which the reference 
range tool we have developed would be particularly useful. It is a 
genetic condition that results in the production of abnormal 
hemoglobin and subsequent chronic hemolytic anemia. As a 
compensatory mechanism for low hematocrit, both high cerebral 
blood velocity and flow have been described.^''"^^ Children with 
SCA are at high risk of both clinical and silent cerebral infarction, 
narrowing of large vessels, and insufficient oxygen/glucose 



delivery, which may result in deficits in brain function even in 
patients without any abnormality seen on r2-weighted MRI. The 
SCA patient included in this study was one such case, in which 
hemoglobin and hematocrit levels were below the normal range, 
but no evidence of infarction was visible on the r2-weighted 
images. However, CBF was well above the age-matched reference 
range in all vascular territories in this patient, being approximately 
double the expected value in a typical 12-year-old patient. This is 
indicative of a compensatory increase in blood flow due to the low 
hematocrit in this patient (it should be noted that quoted values 
of CBF in this patient may be slightly inaccurate, as the literature 
value of Tib used does not account for the low hematocrit in this 
patient). Similarly, the short values of 7 observed in the ACA and 
MCA territories are likely to be due to increased blood flow 
velocities in the feeding arteries, resulting in a reduction in the 
level of bolus dispersion. Reduced values of Ti were observed in 
all vascular territories on the RHS, and bilaterally in the MCA 
territory. This agrees with a previous study, which showed 
reduced values of Ti in the GM of SCA patients, even in the 
scans of patients appearing normal on conventional MRI.'^^ It has 
been suggested that this is indicative of brain damage due to 
moderate hypoxia, and is particularly present in tissues with a 
higher metabolic rate,^^ which could explain the reduction in GM 
Ti values seen in this patient. 



CONCLUSION 

Overall, our results show that during the transition from childhood 
into adulthood, significant changes occur in both and Mq of 
cortical GM and the hemodynamic properties of CBF. Due to the 
strong age dependence of Ti and Mq, if assumed rather than 
measured values are used when calculating CBF from raw ASL 
data (as is often the case in clinical studies, when scan time is 
limited), age-corrected values should be used, as adopting adult 
values of Mq and Ti in equation (1) will lead to a significant 
overestimation of CBF in children. Similarly, once CBF values have 
been calculated, the data should be compared with age-matched 
reference range values, due to the normal variation with age in 
healthy subjects in this age range. In multi-TI ASL studies, normal 
values of T will also vary as a function of age (depending on the 
vascular region), which should be accounted for when designing a 
protocol to investigate pathologies in this region. 
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